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Abstract 


Theoretical  study  of  hypersonic  boundary-layer  receptivity  to  acoustic  waves  scattered  by 
a  wall  waviness  (distributed  roughness)  and  hump  (local  roughness)  is  performed  using  a 
combination  of  asymptotic  and  numerical  techniques.  Analytical  results  have  been 
obtained  for  the  normal-wave  excitation  by  any  combination  of  external  disturbances 
(sound,  vortical  and  entropy  waves)  and  wall  perturbations  (roughness,  suction/blowing, 
temperature  irregularities  etc.).  It  is  shown  that  a  strong  excitation  occurs  in  local  regions 
where  the  flow  disturbance,  which  is  due  to  nonlinear  interaction  of  external  waves  with 
wall-induced  disturbances,  is  in  resonance  with  normal  waves.  Analytical  results  are 
integrated  into  the  computational  module  predicting  the  normal-wave  amplitude 
generated  in  the  resonance  regions.  Calculations  of  the  second-mode  excitation  by 
acoustic  waves  of  different  frequency  and  incidence  angle  were  conducted  for  two- 
dimensional  disturbances  in  the  boundary  layer  on  a  flat  plate  at  Mach  6.  Receptivity 
functions  are  presented  for  the  case  of  both  distributed  and  local  roughness.  The 
theoretical  model  can  be  used  to  predict  initial  amplitudes  of  unstable  normal  waves  as 
well  as  estimate  a  tolerable  roughness  size  associated  with  transition. 
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Nomenclature 


Latin  symbols 

A 

A 


f 

G 

g 

L 

M 


eigenfunction  vector  (direct  problem) 

eigenfunction  vector  (conjugate  problem) 
phase  speed 

frequency  parameter 

streamwise  component  of  mass-flux  disturbance 

receptivity  function 

roughness  shape-function 

streamwise  scale 

Mach  number 


p  =  p*  pressure  disturbance 

p-p*  IP*  -  jM'^p  pressure  disturbance 


q 

X  =  py/y: 

t 

T 

U 

V 

u,v,w 

x,y,z 


X,  =X*/Z‘  =£X 


receptivity  factor 

Reynolds  number 
time 

mean  flow  temperature 
mean  flow  streamwise  velocity 
mean  flow  vertical  velocity 
velocity  disturbance 
Cartesian  coordinates 
slow  variable 


Greek  symbols 

cc,  fi  wavevector  components  in  x-  and  z-direction  respectively 

Y  specific  heats  ratio 
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s  =5* ! C  small  parameter 

A*  boundary-layer  thickness,  J7*(A*)  =  0.99f/* 

s'=yjv:v/u:  boundary  layer  scale 

V  kinematic  viscosity 

p{cc,fi)  F ourier  component  of  the  shape-function  z ') 

CO  angular  frequency 

0  incidence  angle 

9  temperature  disturbance 

’5'  disturbance  vector-function 

F{x,y)  disturbance  amplitude  vector-function 

Subscripts 

0 
a 

ad 
b 

dis 
e 

loc 
w 

Superscripts 

bar 

asterisk 

Abbreviations 

Im  imaginary  part 

Re  real  part 


complex  conjugate  value 
dimensional  value 


resonance  point 
acoustic 
adiabatic  wall 
branch  point 
distributed 

upper  boundary-layer  edge 

local 

wall 
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1 .  Introduction 


This  report  summarizes  the  third  phase  of  our  theoretical  studies  on  hypersonic  boundary- 
layer  receptivity.  The  first  phase  [1]  addressed  receptivity  of  two-dimensiortal  boundary 
layer  on  a  flat  plate  or  sharp  cone  to  entropy/vorticity  disturbances.  It  was  found  that  the 
first  and  second  modes  are  synchronized  with  these  freestream  disturbances  near  lower 
neutral  branch.  Due  to  nonparallel  effects  entropy/vorticity  waves  are  partially  swallowed 
by  the  boundary  layer  and  effectively  generate  the  boundary-layer  instability  in  the 
synchronism  region.  This  mechanism  can  compete  with  the  leading  edge  receptivity  to 
acoustic  disturbances  [2,  3]  in  cases  of  “quiet”  freestream  and  conical  body 
configurations. 

In  the  second  phase  [4],  we  analyzed  receptivity  to  wall-induced  disturbances  such  as 
vibrations,  periodic  suction-blowing  through  holes  and  slots  and  wall-temperature 
perturbations.  This  receptivity  mechanism  is  associated  with  active  (ascend)  flights  of  a 
hypersonic  vehicle,  when  propulsion  system  generates  skin  vibrations.  The  latter  can 
propagate  upstream  and  excite  boundary-layer  instabilities  on  the  forebody  surface. 
Receptivity  to  wall-induced  disturbances  is  also  important  in  developing  active 
techniques  providing  dynamic  control  of  initial  instability  amplitude.  In  these  techniques, 
wall  actuators  generate  disturbances  in  counter-phase  with  “natural”  unstable  waves  in 
order  to  cancel  the  existing  instability.  The  theoretical  analysis  [4]  showed  that 
hypersonic  boundary  layer  is  extremely  receptive  to  wall  vibrations  of  resonance 
firequency  and  wavelength.  It  was  found  that  the  receptivity  function  is  maximal  near  the 
lower  neutral  branch  of  the  second  mode.  This  maximum  is  associated  with  the 
disturbance  spectrum  branch  point  at  which  characteristics  of  the  first  and  second  modes 
are  identical.  Near  this  point,  both  excitation  and  propagation  of  unstable  normal  waves 
are  singular.  With  the  help  of  asymptotic  analysis  we  resolved  singularities  and  evaluated 
the  receptivity  function  in  the  branch-point  vicinity.  This  study  showed  that  a 
straightforward  extrapolation  of  the  receptivity  and  stability  procedures  from  moderate 
supersonic  to  hypersonic  regimes  may  lead  to  wrong  predictions. 
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Another  important  mechanism  of  hypersonic  boundary-layer  receptivity  is  relevant  to  the 
interaction  of  external  unsteady  disturbances  (acoustic,  entropy  and  vorticity  waves)  with 
steady  disturbances  induced  by  spatial  variations  in  surface  boundary  conditions 
(roughness/waviness,  local  heating/cooling,  suction/blowing  etc.).  Early  experimental  [5] 
and  theoretical  [6-8]  studies  of  this  mechanism  have  been  conducted  in  Russia  for 
subsonic  flows.  Their  results  have  been  summarized  in  [9].  Asymptotic  theory  of  the 
Tollmien-Schlichting  wave  excitation  by  acoustic  wave  interacting  with  surface 
irregularities  has  been  developed  in  [10-12]  and  reviewed  in  [12-14],  In  the  last  decade, 
detailed  parametrical  studies  of  this  receptivity  mechanism  have  been  represented  in  [IS¬ 
IS].  Theoretical  predictions  are  in  a  good  agreement  with  the  experiment  [19]  conducted 
on  a  flat  plate  in  low-speed  wind  tunnel. 

Receptivity  of  high-speed  boundary  layers  (in  the  Mach  number  range  from  1.15  to  4.5) 
to  ffeestream  disturbances  scattering  by  the  surface  irregularities  has  been  analyzed  by 
Choudhari  and  Streett  in  [20].  They  considered  generation  of  inflectional  and  higher 
modes  instability  (the  first  and  second  modes)  at  Mach  4.5.  Choudhari  and  Streett  [20] 
reported  that,  although  the  Mach=4.5  boundary  layer  is  nearly  6  times  thicker  than  in  the 
Mach=0  case,  the  instability  amplitudes  generated  by  a  wall  hump  of  fixed  height  are 
only  marginally  smaller  than  those  at  low  speeds.  This  example  indicates  that  receptivity 
of  hypersonic  boundary  layer  to  acoustic  plus  roughness  disturbances  may  play  important 
role  in  the  initial  phase  of  transition. 

In  this  report,  we  perform  theoretical  analysis  of  h5^ersonic  boundary-layer  receptivity  to 
acoustic  waves  scattered  by  surface  waviness  (distributed  exciation)  or  hump  (local 
excitation)  on  a  flat  plate.  In  Section  2,  we  formulate  the  problem  and  show  that  main 
features  of  its  analysis  are  similar  to  the  case  of  receptivity  to  wall  disturbances 
considered  in  [4].  In  Section  3,  we  discuss  the  numerical  results  obtained  for  the  second¬ 
mode  excitation  by  fast  and  slow  acoustic  waves  of  various  incidence  angles  in  the 
boundary  layer  at  Mach  6.  In  Section  4,  we  conclude  the  report  and  discuss  future  effort. 
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2.  Problem  formulation  and  analysis 

In  this  chapter,  we  consider  the  excitation  of  boundary-layer  normal  waves  by  an 
acoustic  wave  which  incidents  on  a  flat  plate  and  interacts  with  its  wavy  surface 
(distributed  receptivity)  or  hump  (local  receptivity).  In  our  analysis,  we  use  the  method 
developed  in  [4]  for  receptivity  to  wall-induced  disturbances. 

2. 1  Distributed  receptivity  on  wavy  surface 

We  consider  two-dimensional  laminar  boundary  layer  of  a  perfect  gas  on  a  hypersonic 
flat  plate.  The  global  streamwise  length  L*  is  assumed  to  be  much  larger  that  the 

boundary  layer  thickness  scale  8*  =  /Ul ;  i.e.  the  ratio  s  =  8*/L*  is  small  and  the 

Reynolds  number  =e"‘  is  large.  The  longitudinal  x*,  normal  y*, 

transversal  z*  coordinates  and  time  t*  are  made  nondimensional  using  5*  and  U]  as 

{x,y,z)  =  {x\  y,  z*)l5\  t  =  tV:/8\ 

Introducing  the  slow  variable  x^  =  x*  I V  =£x  we  can  specify  the  nondimensional  mean- 
flow  velocity  components  {U,  V),  temperature  T,  and  pressure  P  as 


U  =  U* /U*  =U(x„y),  V  sV’/t/;  =£Vo(x,,y),  T  ^T*  IT;  =T(x„y) ,  (2.1) 

P^P*l{p:uf)  =  P{x,). 


The  wall  waviness  is  represented  in  the  form 

yAx,  y,  z)  =  /texp(ia„x  -i-  ifi^z)  +  (c.c) ,  (2.2) 

where  h  =  h*  15*  is  nondimensional  amplitude,  is  wavenumber  vector.  It 

is  assumed  that  the  wavelength,  X„  =  2kIk^=  0(1) ,  is  of  the  order  of  boundary  layer 
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thickness  A’  (f/*  =0.99U*  at  3/*  =  A*).  The  waviness  amplitude  h*  is  smaller  than  the 

viscous  sublayer  thickness;  i.e.  h*  <  .  In  this  case,  the  wavy  wall  generates  small 

steady  disturbances,  which  are  governed  by  linearized  Navier-Stokes  equations  [21]. 

An  external  acoustic  wave  of  prescribed  angular  frequency  co  =  co*5*  I U*  and  wave- 
numbers  ag=alUl/5\  incidents  upon  the  plate  as  schematically 

shown  in  Fig.  1 .  Its  physical  quantities  can  be  represented  in  a  traveling-wave  form.  For 
example,  the  incident  wave  pressure  is  expressed  as 

P  =  Pa.mc  exp(ia^x  -1-  iP^z  +  y-i(Ot)  +  (c.c) .  (2.3) 


The  vertical  wave-number  component,  ,  is  determined  from  the  acoustic  dispersion 

relation  (0,  M) ,  where  M  is  freestream  Mach  number.  The  incident 

wave  is  refracted  by  the  boundary  layer,  reflected  by  the  wall  and  scattered  by  the  wall 
waviness.  Interacting  with  the  waviness-induced  disturbance  the  acoustic  wave  generates 
boundary-layer  normal  waves  including  the  second  mode.  Our  objective  is  to  determine 
the  receptivity  function  coupling  characteristics  of  the  incident  acoustic  wave  and  the 
wall  waviness  with  the  normal-wave  amplitude. 


Using  the  approach  of  [4]  we  represent  the  unsteady  flow  field  as  the  vector  function 


,7,.  ,  .  du  „  -  d9  dw.T 

^'(x,y,t)  =  (u,—,v,p,e,—,w,—)  , 
oy  dy  ay 


(2.4) 


where  u,  v  and  w  are  velocity  components;  p  =  p  /(p^U^')  is  pressure  and  6  is 
temperature.  In  the  following  analysis  we  will  also  use  the  nondimensional  pressure 
p  =  p*  /  P*  =yM'p .  The  vector-function  (2.4)  is  expressed  in  the  form 
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(2.5) 


't'  =  y)  +  hP^ix^ ,  y)exp(fa^j:  +  ifi^z)  + 

+  (Xi ,  y)cxp(ia„x  +  ip^z  -  ica)  +  F(x,  y,  z)  exp(-i(ot)  +  (c.c) 

Here  the  first  term  describes  the  boundary-layer  mean  flow  on  smooth  surface;  the 
second  term  corresponds  to  the  stationary  disturbance  induced  by  the  wall  waviness.  The 
third  term  corresponds  to  the  acoustic  disturbance,  which  is  a  combination  of  incident 
and  reflected  waves  plus  viscous  components  (vorticity  and  entropy  waves)  associated 
with  the  interaction  of  the  acoustic  component  with  the  boundary-layer  flow.  Definition 
of  the  amplitude  depends  on  the  acoustic  wave  normalization.  Hereafter  we  normalize 

the  vector-fiinction  F^(y)  by  the  condition:  (pressure  amplitude  p  of  the  incident 
acoustic  wave)=l  at  the  upper  boundary-layer  edge;  then  =  p*a  j^l  P* . 

The  fourth  term  in  Eq.  (2.5)  represent  disturbances  generated  in  the  boundary  layer  due 
to  nonlinear  interaction  between  the  acoustic  wave  and  the  waviness-induced 
perturbation.  It  is  assumed  that  disturbance  amplitudes  are  small,  « 1  and  ^  « 1 , 

and  the  higher  order  terms  such  as  e^h ,  are  neglected;  i.e.  the  nonlinear  interaction 
is  considered  to  be  weak. 


Substituting  (2.5)  into  Navier-Stokes  equations  we  obtain  the  following  problems,  which 
can  be  written  in  the  matrix-operator  form. 

For  stationary  disturbances  induced  by  wall  waviness: 


H{y,d  ,x^,(O  =  0,a^,PjF^=0, 


(2.6a) 


KM„0)  =  -—{x,Si),  F^^{x„Q>)  =  0,  F^,{x,,0)  =  -—{x,fi),  F„7(Xi,0)  =  0,  (2.6b) 


<  <»  ,  y  — >  oo . 


(2.6c) 
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Boundary  conditions  (2.6b)  have  been  obtained  by  the  Taylor  series  expansion  of  the  no¬ 
slip  conditions,  m  =  v  =  w  =  0 ,  and  the  wall-temperature  condition,  6  =  T„,z.i  y  =  y,, .  In 

a  similar  way  we  can  consider  a  wide  class  of  inhomogeneous  boundary  conditions  on 
the  wall  surface  such  as  suction/blowing  and  heating/cooling.  These  variations  do  not 
affect  the  analysis. 

For  acoustic  disturbances: 

H(y,dy,x^,0},a^,p^)F^  =0, 

F^;(Xi,0)  =  F„3(Xi,0)  =  F^jCxpO)  =  F„^(x,,0)  =  0, 

(pressure  amplitude  of  incident- wave  component)=l,  y 

For  boundary-layer  disturbances: 

(y,  ^3,,  Xp  ,  z,  dj ,  ty)F 

=  hs^Gi y, Xj ,  F^ ; 0},a„ exp[i(a^  +(xjx  +  i(p„  +  )z] ,  (2.8a) 

Fi(x„0)  =  F3(XpO)  =  F5(XpO)  =  F,(XpO)  =  0,  (2.8b) 

Fj  <  00  ,  y  — »  00  .  (2.8c) 

The  right-hand  side  of  Eq.  (2.8a)  describes  nonlinear  (quadratic)  interaction  of  the  wall- 
induced  steady  disturbance  with  the  acoustic  wave. 

In  Eqs.  (2.6a)  and  (2.7a),  the  operator  H  can  be  expressed  in  the  standard  form 

H=-^-H,(y,x,-,R,Q),a,p),  (2.9) 

dy 

where  the  matrix  has  dimension  8  x  8 .  Its  elements  depend  on  the  mean-flow  profiles 

(2.1),  disturbance  parameters  a,  P,  O)  and  Reynolds  number  R  =  yju*x*  /  v*  .  Explicit 
form  of  the  matrix  is  given  in  Appendix. 


(2.7a) 

(2.7b) 

(2.7c) 
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We  consider  a  partial  solution  of  the  problem  (2.8a)-(2.8c)  as  a  decomposition  of  the 
discrete  modes  with  eigenvalues  ai^(x^,P,co) 

^  =  he  [^*0  (^1  ’  y)  +  (^1 ,  y)  +  •  •  •]  exp(/e  +  ifiz) , 

k 

5,  =  ^a,{x„(o,^)dx, . 

■*-1j 

If  the  nonlinear  term  of  Eq.  (2.8a)  is  not  in  resonance  with  the  normal  waves;  i.e. 

and/or  +  then  the  normal- wave  excitation  is 

exponentially  small  and  can  be  neglected.  In  this  case,  we  obtain  the  standard  eigenvalue 
problem 


{S  „  1 

j 

^ko  ®  > 

(2.10a) 

(^*0)13.5,7  — 0, 

(2.10b) 

If^o  ^Oaty-^00. 

(2.10c) 

Its  solution  describes  propagation  of  normal  waves  in  the  boundary  layer  far  from 
resonance  regions  and  can  be  expressed  as 

Pico=c,(x^)A,(x^,y,a,),  (2.11) 

where  are  eigenfunctions  normalized  by  a  certain  condition.  For  example,  the 
pressure  disturbance  amplitude  is  constant  on  the  wall:  A^4(x^,0,a^)  =  1 .  The  amplitude 
functions  q  (jc, )  are  unknown  on  this  stage. 

If  the  nonlinear  term  is  in  resonance  with  one  of  the  normal  waves  (say  the  wave  with 
eigenvalue  )  at  some  point  x^  =  x,o , 


(2.10a) 

(2.10b) 
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a„ I3)  =  a„+a^  and  P  =  +  13^, 


(2.12) 


then  this  normal  wave  is  effectively  generated  by  the  acoustic  wave  due  to  relatively 
strong  sound-waviness  interaction  near  the  resonance  point  .  We  consider  the  case 

when  the  normal  wave  is  unstable  and  the  resonance  point  is  close  to  the  lower  neutral 
branch.  Scheme  of  the  resonant  excitation  region  is  shown  in  Fig.  2.  Analysis  of 
governing  terms  in  Eq.  (2.8a)  indicates  that  the  resonance  excitation  occurs  in  the 
relatively  small  region  |j:j  -Xjo|  =  ,  where  we  can  introduce  the  inner  variable 


(2.13) 


We  assume  that 

P  =  Pw+Pa  cc^+oCa=(Xo+  ye"- ,  (2.14) 

where  =a„{x^Q,(0,P) ,  and  the  real  parameter  y  characterizes  resonance  detuning. 
Now  the  inner  solution  of  the  problem  (2.8a)-(2.8c)  can  be  expressed  in  the  form 

F  =  +  A(l.  3^)  •  •  .]exp(/£-‘^^a<,|  +  i^z) ,  (2. 15) 

where  Afy)  =  A„(ji:jg,y,ao)  is  normal- wave  eigenfunction  at  the  resonance  point. 
Expanding  the  operator  in  the  vicinity  of  a:  =  Xjg  we  obtain  the  following  problem 
for  the  vector-function 


H  q(^X'^q,(Xq) 

\dy 


\ 

A 


-i  ^^10’  -^(^io>  «o)A  +  GPa  exp(/)^) , 

(2.16a) 

(A)i,3.5.7=0  aty  =  0,  (2.16b) 
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— >  0  at  j  -»  00 . 


(2.16c) 


This  problem  is  similar  to  the  problem  (2.20)  in  [4],  It  has  non-trivial  solution,  if  its  right- 
hand  side  is  orthogonal  to  the  eigenfunction  .B,,  =  jB„(x,p,j,ao)  of  the  conjugate 
problem.  This  condition  leads  to  the  following  equation  for  the  amplitude  coeflScient  c(^) 


^-/6^c  =  ^exp(/>^), 


,  da„  ,  ^ 


q  =  i- 


5  ^0  /  \ 

^  a  (^10>^o)^0  ^ 


where  the  scalar  product  is  defined  as 


< 


0  J=^ 


(2.17a) 

(2.17b) 

(2.17c) 


(2.18) 


the  upper  bar  denotes  complex  conjugate  values.  Nonzero  elements  of  the  nonlinear 
interaction  matrix  G  are  represented  in  Apendix.  If  the  normal  wave  has  zero  amplitude 
upstream  from  the  resonance  region,  then  the  boundary  condition  for  Eq.  (2.17a)  is 
specified  as 


c(^)  0  at  ^  -00 .  (2. 19) 

Problem  (2.17a)-(2.19)  coincides  with  the  problem  (2.21a)-(2.22)  of  [4],  if  we  replace  the 
numerator  (.Bo,Oo)  in  Eq.  (2.21c)  by  the  scalar  product  < Bo,G(x^q)F^  >.  Using  this 
similarity  we  obtain  the  normal-wave  amplitude  downstream  from  the  resonance  region 
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2b 


J 


IS 


•'i 

ja„(x,)c&, 


*10 


for  X,  >  x„ 


(2.20) 


This  equation  shows  that  the  normal-wave  amplitude  is  proportional  to  (Z*  / , 
whereas  the  resonant  region  length  is  /*  =  .  As  the  Reynolds  number  increases, 

nonparallel  effects  (which  disturb  the  resonance  condition)  decrease  and  the  relative  length 
rtS”  increases.  In  the  parallel  mean-flow  limit,  oo ,  the  resonance  region  is 
infinitely  large  and  the  normal  wave  amplitude  tends  to  infinity.  In  this  case,  the  upper 
limit  of  the  disturbance  amplitude  is  determined  by  nonlinear  effects  of  the  higher  order 
approximation.  These  conclusions  are  similar  to  those  formulated  in  [4]  for  receptivity  to 
wall  disturbances. 


Note  that  nonparallel  effects  make  the  receptivity  mechanism  to  be  local  in  units  of  Z* 
even  for  distributed  roughness.  However  it  can  be  treated  as  a  distributed  mechanism  in 
units  of  <5* .  Our  choice  of  the  term  “distributed  receptivity”  is  a  matter  of  convention. 

2.2  Local  receptivity  on  hump 

We  consider  a  spanwise  row  of  local  humps  on  a  flat  plate  surface.  For  brevity  we  analyze 
one  spanwise  harmonic  of  this  periodic  structure,  and  specify  its  shape  as 

yAx,y,z)  =  hg(x')exp(ifi„z)  +  (c.c),  x'  =  x-x^  =£-'(x,  -Xi^).  (2.21) 


The  shape-function  is  represented  by  the  Fourier  integral 


_i_ 

Itt 


-l-co 


j  p(ajexp(ia„x')da 


w  ? 


(2.22) 


which  can  be  expressed  in  the  form 
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2  +00 

S  =  J ^ '"V)  exp(ir^)dr .  (2.23) 


Here  cc„o  =  a„(x^Q,co,P)-a^  is  wavenumber  component  fitting  to  the  resonance 
condition  (2.12),  and  y  characterizes  resonance  detuning  in  accordance  with  Eq.  (2.14). 

Using  the  relation  (2.20)  for  the  normal-wave  amplitude  excited  by  one  harmonic  and 
integrating  over  y  we  obtain 


Pn  ='^AP(«woK«3^)exp 


X 

ija„(x')dx'  +  ifiz 


\  0 


(2.24) 


This  expression  is  similar  to  Eq.  (2.33)  of  [4]  for  a  local  wall  forcing.  Important  that  the 
function  ^(x,o,aQ,Qr^,^^)  does  not  depend  on  the  hump  shape  g(x') .  If  this  function  is 

calculated  once  and  for  all,  then  we  can  evaluate  the  normal  wave  amplitude  generated  by 
the  sound  wave  interacting  with  a  wall  hump  of  arbitrary  shape  (within  the  limit  of  our 
assumptions). 


The  acoustic  wave  interacting  with  a  local  hump  generates  the  boundary-layer  normal 
wave  of  a  finite  amplitude  even  in  the  parallel  mean  flow.  This  is  due  to  the  fact  that  only  a 
small  portion  of  the  hump  spectrum,  (a^  -a^o)  =  is  involved  into  the  resonance 

mechanism.  As  ^  ^  0,  the  resonant  excitation  increases  proportionally  to  ,  whereas 
the  active  spectrum  region  is  narrowed  down  proportionally  to  .  Because  both  trends 
compensate  each  other,  the  normal-wave  amplitude  tends  to  a  finite  value. 

Note  that  Eqs.  (2.20)  and  (2.24)  hold  for  any  combinations  of  fireestream  disturbances  and 
wall  perturbations.  For  example,  if  we  replace  the  acoustic  wave  by  vortical  (entropy) 
wave,  then  these  expressions  provide  the  normal  wave  amplitude  generated  by  vortical 
(entropy)  wave  interacting  with  a  surface  irregularity.  The  normal  wave  can  also  be  of 


21 


different  nature  such  as  Tollmien-Schlichting  wave,  second  mode,  cross-flow  vortices, 
Gortler  vortices  etc.  All  these  cases  can  be  treated  using  the  same  computational  module, 
which 

1 .  solves  Problem  (2.6a)-(2.6c)  for  wall-induced  disturbance; 

2.  solves  Problem  (2.7a)-(2.7c)  for  external  (freestream)  disturbance; 

3.  solves  the  eigenvalue  problem  (2.10a)-(2. 10c)  and  its  conjugate  problem  to  obtain  the 
normal-wave  eigenfunctions  A^(x^,y)  and  B^(x^,y); 

4.  evaluates  the  scalar  products  of  Eqs.  (2. 17c); 

5.  predicts  the  normal-wave  amplitude  using  Eq.  (2.20)  for  distributed  receptivity  or  Eq. 
(2.24)  for  local  receptivity. 

Such  a  module  was  developed  and  coupled  with  our  stability  solver. 

2.3  Receptivity  near  spectrum  branch  points 

Foregoing  analysis  is  valid  for  “simple”  discrete  spectrums;  i.e.  eigenvalues  of  the 
boundary-layer  modes  are  assumed  to  be  different.  This  is  typical  for  subsonic  and 
moderate  supersonic  boundary  layers  as  well  as  hypersonic  boundary  layers  on  adiabatic 
walls.  The  analysis  of  [1,  2,  4]  showed  that  the  first  and  second  modes  have  branch  points 
which  are  close  to  the  real  parameter  space  in  the  case  of  low  wall  temperature  ratios.  The 
upstream  branch  point,  x,  =  Xjj ,  is  located  near  the  lower  neutral  branch  of  the  second 
mode  whereas  the  downstream  branch  point  is  close  to  the  upper  neutral  branch.  In  the 
branch-point  vicinity  (say  for  the  upstream  point  ),  the  eigenvalues  behave  as 

cc^^2  —  ^  Xj  ,  (^*^^) 

where  =  a,(x,6,<»)  =  a^ix^i,^). 
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The  following  analysis  is  focused  on  the  upstream  branch  point.  The  receptivity  factor  q 
in  Eq.  (2.17c)  tends  to  infinity  as  the  resonance  point  tends  to  the  branch  point, 
■*10  ^  -*16  •  This  singularity  is  due  to  the  fact  that  the  scalar  product 

^0,as  (2.26) 

The  eigenvalue  derivative,  which  is  involved  into  Eq.  (2.20)  for  distributed  receptivity, 
behaves  as 


^  da 


(•*10 A)’  ^0 


b  =  -^(^lo)  ^  ,  as  ^  .  (2.27) 

This  singularity  has  been  studied  in  [4]  for  the  case  of  receptivity  to  wall-induced 
perturbations.  Since  the  analysis  of  [4]  is  easily  extended  to  the  problem  considered  in 
this  report,  we  will  briefly  indicate  key  points  of  this  extension  and  summarize  the 
results.  For  brevity,  we  consider  two-dimensional  normal  waves  with  P  =  0 . 

Near  the  branch  point.  Mode  1  (with  the  eigenvalue  a^(x^,o)))  strongly  interacts  with 

Mode  2  (with  the  eigenvalue  a,(x,,fi)))  due  to  nonparallel  effects.  This  interaction  can 

be  described  by  the  two-mode  approximation;  Le.  the  disturbance  amplitude  from  Eq. 
(2.5)  is  expressed  in  the  form 

F(Xi,y)  =  c,(Xi)/4i(x,,y)exp(e“‘r5,)-i-c,(Xj)4(x,,y)exp(£“‘/52) .  (2.28) 

Far  from  the  resonance  regions,  the  amplitude  coefficients  Ci  2(a:,)  are  solutions  of  the 
ODE  problem 


dq 

dx. 


I  _ 


=  +  CjW,,  exp[e  '/(5,  -  Sj)] , 


(2.29a) 
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(2,29b) 


^  +  qW,,  exp[e-'/(5,  -  5, )] , 

dx^  '  ~  ~ 

~  ^2s  ’  (2.29c) 

where  and  c,^  are  prescribed  values  at  a  certain  initial  point  ;  the  matrix  elements 
are  defined  as 


W  =■ 


d<X  dXi 


SA, 

da 


7  =  1,2,  k  =  \,2; 


(2.30) 


where  the  matrix  i/,  is  relevant  to  nonparallel  effects.  Its  elements  depend  on  the  mean- 
flow  profile  derivatives  dU  /  dxj ,  dT  /  dxj  and  vertical  velocity  .  Near  the  branch  point, 
the  matrix  elements  (2.30)  behave  as 


(_1);>C-1 

^jk  - TT’  ■ 


4(ar, -Xj^) 


(2.31) 


That  leads  to  the  following  singularity  of  the  amplitude  coefficients 


c 

'^1,2 


(x,-xj 


1/4 


-I-...,  as  JCj 


(2.32) 


where  constants  C,  ,  depend  on  the  normalization  of  eigenfunctions  Aj  2(.Xi,y). 

In  the  region  x^  -x^^  =0(e^'^),  the  mode  decomposition  (2.28)  is  replaced  by  the  local 
expansion  with  the  inner  variable  ^  =  e~''^(x^  -  x^^) 
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F  =  £~''^[co(0\(y)  +  £~''"coi(0\i(y)  +  ■■■]e^p(i£~'S,  +  ■  (2-33) 


Here  Cq(^)  is  solution  of  Airy  equation 


dC- 


-A^Cco=0. 


(2.34) 


Matching  the  inner  and  outer  solutions  we  can  establish  an  exchange  rule  in  the  branch¬ 
point  vicinity  [2,  4],  The  local  expansion  (2.33)  indicates  that  the  disturbance  amplitude 
has  a  peak  in  the  branch-point  influence  domain  ^  =  0(1) ,  with  maximum  amplitude 

being  of  the  order  of  e  .  To  separate  this  effect  from  the  receptivity  mechanism  we 
need  to  establish  a  relationship  between  the  constants  (instead  of  the  amplitude 

functions  Cj , )  and  characteristics  of  incident  acoustic  wave  plus  wall  disturbance. 

Wavy  wall 

If  the  incident  acoustic  wave  interacts  with  the  wall  waviness,  then  the  unstable  normal- 
wave  amplitude  is  expressed  as 


Fn 


-1/2 


^10 


,  (2.35) 


•*■1  ^  •'‘10  • 


Whereas  the  local  amplitude  coefficient  of  the  unstable  mode  n  is 


/  d//  ' 

(^0’ 


as  Xio  ->  Xj* . 


(2.36) 
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Equation  (2.36)  shows  that  distributed  receptivity  has  no  singularity  at  the  spectrum 
branch  point.  In  accordance  with  Eq.  (2.32),  the  disturbance  peak  is  totally  due  to  a 
singular  behavior  of  the  unstable  mode  itself.  If  one  measures  the  disturbance  amplitude 
Q„  at  a  certain  fixed  point  downstream  from  the  branch  point  for  various  locations  of  the 

resonant  point  then  he  will  not  observe  significant  changes  in  the  distribution 
Q„(x^o)  near  the  point  Xjq  =  .  This  conclusion  is  similar  to  that  formulated  in  [4]  for 

receptivity  to  wall-induced  disturbances. 

Local  hump 

If  the  hump  is  far  from  the  branch  point,  then  the  normal-wave  amplitude  is  determined 
from  Eq.  (2.24).  As  the  hump  approaches  the  branch  point,  the  receptivity  factor  behaves 
as 


q-^ 


,  as  XiQ  ->  x,. . 


(2.37) 


Then  the  amplitude  coefficient  from  Eq.  (2.32)  is  expressed  in  the  branch  point  vicinity  in 
the  form 


U  BH 


Ap(ao) 


-  \  2(Xjq  Xjj) 


1/4 


,  as  Xjo  — >  X] 


'•1*  ■ 


(2.38) 


Local  receptivity  increases  proportionally  to  (Xjo-x,^)  as  the  hump  approaches  the 

branch  point.  This  effect  is  due  to  the  increase  of  the  active  wavenumber  range  of  the 
hump  spectrum.  If  the  hump  is  far  from  the  branch  point,  then  this  range  is  estimated  as 
(a^+a^  -Uq)  =  0{e^'~) .  If  the  hump  is  placed  in  the  region  (Xjo  -x,j,)  =  Oie''^) ,  then 

the  active  wavenumber  range  becomes  of  the  order  of  -i-a^  -cTo)  =  In  this 

case,  the  amplitude  coefficient  is  estimated  as  =  0(e^/i£“^'®) . 
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Analysis  of  the  inner  solution  (2.33)  and  its  matching  with  the  outer  solution  lead  to  the 
following  expression  for  the  normal  wave  amplitude 


Ap(a,)r(l/3) 


f  Q  ^  5 


2V7r(Xi-XiJ 


1/4 


V  ^  / 


4„(Xi,y)exp 


^10 


(2.39) 


where  « Xj -a:,^  « 1;  subscript  “/>”  denotes  quantities  evaluated  at 

Gamma-function  r(l  / 3)  =  2.678938....  From  this  expression  and  Eq.  (2.32)  we  obtain 
the  amplitude  coefficient 


Cn  =£ah\ 


1/6 


Ap(«,)r(l/3)<g„G(A,JF,> 


B 

\ 


(2.40) 


Summarizing  we  conclude  that  local  receptivity  strongly  depends  on  the  distance 
between  the  hump  and  branch  point.  As  Aio->Ai*,  the  amplitude  coefficient  C„ 

increases  proportionally  to  |a:io-a,^|  and  attains  its  maximum  value  of  the  order  of 

s^he  at  the  branch  point  a^j  =  Aj^,  .  If  one  measures  the  disturbance  amplitude  at  a 
certain  fixed  point  downstream  from  the  branch  point  for  various  hump  locations,  then  he 
will  observe  a  local  peak  in  the  distribution  Q^iXy^)  near  the  point  Aj^  =  a,*  .  Because  the 

branch  point  is  close  to  the  lower  neutral  branch,  this  peak  will  be  enhanced  by  the 
exponential  growth  of  unstable  disturbance.  This  result  is  consistent  with  that  obtained  in 
[4]  for  receptivity  to  wall-induced  disturbances. 
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3.  Numerical  results 

Receptivity  calculations  were  conducted  for  the  boundary  layer  on  a  flat  plate  at  local 
Mach  number  M  =  6,  stagnation  temperature  =  390°  F,  and  wall  temperature 

T*  =80°  F  that  corresponds  to  the  wall  temperature  ratio  T*/Tq  =0.6356.  Dr.  Roger 

Kimmel  indicated  this  case  as  a  good  candidate  for  the  Mach  6  wind  tunnel  experiment. 
It  is  assumed  that  the  fluid  is  a  perfect  gas  with  a  constant  specific  heats  ratio  y  =  1.4 
and  Prandtl  number  Pr  =  0.72 .  The  viscosity-temperature  dependency  is  approximated 
by  Sutherland’s  law 

(r-HS) 

where  5  =  110/7’/  for  air  temperature  measured  in  degrees  Kelvin;  )U  =  jU*/jU*  is 
nondimensional  viscosity.  Hereafter  the  global  length  scale  is  L’  =x*,  the  local  scale  is 
self-similar  boundary-layer  scale  S*  =-yJv*x*  /U* ,  the  Reynolds  number  is  =  R  = 

yjuy/v*  and  £=S* /V  =R-\ 


3.1  Mean  flow 

We  neglect  effects  of  hypersonic  viscous-inviscid  interaction  and  assume  that  flow 
characteristics  at  the  upper  boundary-layer  edge  coincide  with  the  ffeestream 
characteristics.  The  laminar  boundary-layer  flow  is  approximated  by  the  self-similar 

profiles  U{y)  and  7’(y),  with  the  boundary-layer  variable  y  =  y* /-^v^x*  lU*  .  These 
profiles  and  their  first  derivatives  are  shown  in  Figs.  3  and  4.  In  this  case,  the  boundary- 
layer  thickness  is  A*  =l6J3-y]vlx*  /U*  (here  U\^)  =  0.99U*),  and  the  displacement 

thickness  A*^  =  13.64-y/v*x*  /t/J .  Figure  4  shows  that  the  wall  temperature  is  slightly 
below  the  adiabatic-wall  temperature. 
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3.2  Second  mode 

Characteristics  of  two-dimensional  normal  waves  have  been  calculated  for  the  second- 
viscosity  coefficient  Jc  =  1.2  ( A:  =  0  corresponds  to  Stokes).  In  Figs.  5  and  6,  the 
eigenvalue  distributions  are  shown  for  various  frequency  parameters 

F  =  (o*v*/uf.  Figures  7  and  8  illustrate  the  phase-speed  distributions  c,(/?)  = 
Re(6)/ a,)  and  the  branch-point  locations  respectively.  The  second-mode  phase  speed  is 
smaller  than  the  freestream  speed  and  larger  than  the  phase  speed  =1-1/ M  of  the 
upper  limit  for  slow  acoustic  waves.  Local  maximums  (minimums)  of  c^iy)  correlate 
with  the  upstream  (downstream)  branch  point  loci.  Note  that  the  branch  points  are  close 
to  neutral  branches  of  the  second  mode.  In  accordance  with  the  analysis  [1,  2],  this 
spectrum  topology  is  typical  for  hypersonic  boundary  layers. 

Figure  8  shows  that  the  upstream  branch  points  are  below  the  real  axis  of  complex  R , 
and  the  downstream  branch  points  are  above  this  axis.  In  what  follows,  we  assume  that 
the  branch-point  influence  domains  are  far  enough  from  the  real  axis  and  use  the 
theoretical  model  of  Section  2  relevant  to  a  simple  discrete  spectrum. 

Distributions  of  the  pressure  amplitude,  p,(y)  =  )M"|A24(y)| ,  and  streamwise  mass-flux 
amplitude. 


/2(y)  = 


T 


■h 


(3.2) 


are  shown  in  Fig.  9  for  the  point  F  =  8.xl0^,  R  =  1000  located  near  the  lower  neutral 
branch.  The  eigenfunction  is  normalized  as  /),(0)  =  1.  The  streamwise  mass-flux 

amplitude  has  a  typical  sharp  maximum  in  the  critical  layer,  where  the  disturbance  phase 
speed  is  close  to  the  mean-flow  velocity. 
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3.3  Waviness-induced  perturbations 

The  pressure  contours  Re[p„(y,a^)]  =  Re[F„^{y,a^)]  =  const  are  shown  in  Fig. 

10  for  the  Reynolds  number  /?  =  1000.  It  is  seen  that  the  wavy  wall  generates  steady 
Mach  waves  outside  the  boundary  layer.  Disturbances  of  maximum  amplitude  are 
observed  at  oc^  -0.75.  In  the  long-wave  ->0)  and  sort-wave  (a„  limits, 

pressure  amplitude  vanishes.  A  major  portion  of  the  boundary  layer  is  weakly  perturbed. 
The  first  maximum  of  pressure  fluctuations  is  located  near  the  upper  boundary-layer 
edge,  where  the  second-mode  critical  layer  is  normally  observed.  Distributions  of 
Re[p„(y)]  are  shown  in  Fig.  1 1  for  =  0.001, 0.05  and  0.1.  It  is  seen  that  the  pressure 

disturbance  amplifies  passing  through  the  laminar  boundary  layer.  Their  amplitude  in 
external  inviscid  flow  is  essentially  higher  than  that  within  the  boundary  layer. 

Contours  of  the  streamwise  mass  flux  component. 


Rel/vvCy^a)]  =  Re 


W'Ka-' 


w5 


1 

—  =  const , 
T 


(3.3) 


are  shown  in  Fig.  12  for  the  Reynolds  number  R  =  1000.  This  pattern  is  qualitatively 
similar  to  the  pressure  pattern  in  Fig.  10  besides  a  sharp  peak  (dark  blue  region)  at  small 
wavenumbers.  This  peak  is  close  to  the  upper  boundary-layer  edge. 

Distributions  Re[/„(y,a)]  at  the  wavenumbers  a  =  0.001,  0.05  and  0.1  are  shown  in 

Fig.  13.  Similar  to  the  case  of  pressure  disturbances,  the  mass-flux  amplitude  outside  the 
boundary  layer  is  essentially  higher  than  that  observed  within  the  boundary  layer. 

Contours  of  the  static  pressure  amplitude,  |p„(y,a)|  =  const,  and  the  streamwise  mass- 
flux  amplitude,  |/„(y,ct)|  =  const,  are  represented  in  Fig.  14  and  15  respectively.  These 
data  give  another  illustration  of  mentioned  above  features  of  the  waviness-induced 
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disturbance.  In  addition,  it  is  seen  that  both  distributions  have  a  flat  maximum  near 
a  »  0.075. 

3.4  Acoustic  waves 

Schematic  pattern  of  acoustic  waves  generated  by  a  cylindrical  source  in  supersonic  flow 
is  shown  in  Fig.  16.  The  incident  wave  front  is  tangent  to  the  flat  surface  at  Point  B.  Near 
this  point,  acoustic  waves  propagate  upstream  and  downstream  with  very  high  phase 
velocity  and  have  the  incidence  angle  0  =  90°.  In  Region  AB,  the  phase  speed  of 
upstream  front  rapidly  increases  from  -  «>  (at  Point  B)  to  0  (at  Point  A,  where  the  Mach 
wave  crosses  the  plate  surface).  In  this  relatively  short  region,  acoustic  waves  propagate 
upstream.  At  the  turning  point  A,  wave  fronts  stagnate  and  form  an  oscillating  Mach 
wave.  Then  they  turn  back  and  move  downsteam  as  slow  acoustic  waves  with  the  phase 
speed  0<c„<  1-1/M.  The  downstream  fronts  propagate  with  the  phase  speed 

decreasing  from  +  «>  (at  Point  B)  to  =  1  + 1  /  M  (in  far-downstream  region).  These 
fronts  form  fast  acoustic  waves. 

To  identify  incident  and  reflected  waves  we  consider  the  acoustic  dispersion  relation  in 
the  inviscid  approximation.  The  pressure  disturbance  can  be  expressed  as 

Pa(x,  y)  =  expiia^x  +  ik^y-  icat)  +  p_  &xp{ia^x-ik^y -i(ot) ,  (3.4) 

where 

K  =  ylMHoc^-co)-  -a;  ,  Re(k^ )  > 0 .  (3.5) 

Analysis  of  the  acoustic-front  kinematics  shows  that 

•  for  counter-flow  waves  (propagating  upstream  in  Region  AB)  and  fast  waves,  p_ 
corresponds  to  incident  wave  and  p^  to  reflected  wave; 

•  for  slow  waves,  p^  corresponds  to  incident  wave  and  p_  to  reflected  wave. 
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This  switch  from  one  component  to  another  is  due  to  supersonic  convection  of  acoustic 
fronts  as  schematically  shown  in  Fig.  16.  It  is  convenient  to  express  the  acoustic 
wavenumber  as  a  function  of  the  incidence  angle  0.  From  the  dispersion  relation  (3.5) 
we  obtain: 

For  slow  waves 


^  _  mMcosQ  ^  m  1 

Mcos0-l’  a„  Mcos0 


,  0  <  0  <  0, 


(3.6) 


where  0^,  =  arccos(l  / M)  corresponds  to  the  Mach  wave  incident  angle. 


For  fast  waves 


G)Mcos0  ,  1  ^  ^ 

cc„  =— - 1 — =1  +  tt - IT’  0<Q<k /2 


MCOS0  +  1 


Mcos0 


(3.7) 


For  counter-flow  waves 


0)M  COS0 


M  COS0- 


1 


Mcos0 


,  0„  <Q<kI1. 


(3.8) 


In  all  figures,  presented  hereafter,  the  acoustic  wave  incidence  angle  0  is  measured  in 
degrees. 

At  h)^ersonic  speeds.  Region  AB  with  counter-flow  waves  is  very  short  and  can  be 
treated  as  a  local  region  concentrated  near  the  line  of  intersection  between  the  Mach 
wave  and  the  body  surface.  In  this  report,  we  do  not  consider  receptivity  to  disturbances 
of  this  type.  We  will  focus  on  the  case  when  the  plate  is  inside  the  Mach  cone  and  its 
surface  is  radiated  by  slow  and/or  fast  acoustic  waves. 
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A  good  guidance  to  the  interaction  of  acoustic  waves  with  a  hypersonic  boundary  layer  is 
provided  by  the  short-wave  (WKB)  analysis  [22,  23].  Schematic  picture  of  this 
interaction  is  shown  in  Fig.  17  for  slow  waves  with  the  phase  speed  a(0)  <c^  <U^-a^, 

where  a(y)  is  local  sound  speed.  In  the  lower  wave-guide,  acoustic  rays  are  reflected  by 
the  wall  and  refracted  near  the  lower  sonic  line,  U (y)  =  -  a(y) .  In  the  upper  half- 

space,  incident  acoustic  waves  are  refracted  near  the  upper  sonic  line,  U(y)  =  c„+a(y). 

In  the  quiet  zone  between  these  two  lines,  the  disturbance  amplitude  decays 
exponentially  with  the  distance  from  the  sonic  lines.  These  qualitative  features  of  the 
acoustic  field-boundary  layer  interaction  are  used  in  the  following  discussion  of 
numerical  data. 

Slow  waves 

Calculations  have  been  performed  for  slow  acoustic  waves  with  the  dispersion  relation 
(3.6)  at  frequency  F  =  8.x  10”^  and  Reynolds  number  R  — 1000 .  The  critical  incidence 
angle  is  0^^  =  80.4°  at  the  Mach  number  M  =  6 .  Figure  1 8  shows  contours  of  the 
pressure  amplitude,  \p^(y,&)\  =  const .  Figure  19  represents  the  function  |/7^(y)|  at 

various  incident  angles.  Pressure  waves  weakly  penetrate  into  the  boundary  layer  because 
its  major  portion  is  covered  by  the  quiet  zone  (see  Fig.  17).  Due  to  interference  between 
reflected  and  incident  waves,  pressure  amplitude  oscillates  from  0  to  approximately  2 
outside  the  boundary  layer;  i.e.  the  reflection  coefficient  is  close  to  one  and  the  boundary 
layer  is  passive  with  respect  to  pressure  waves.  A  flat  local  maximum  of  the  pressure 
amplitude  is  observed  in  the  boundary  layer  at  the  incident  angle  0  ®  40° .  As  0  ->  0„, 
the  disturbance  wavelength  decreases,  the  phase  speed  ^  0,  and  the  upper  sonic  line 

moves  toward  the  wall  allowing  pressure  waves  penetrates  deeper  into  the  boundary 
layer. 

The  streamwise  mass-flux  contours  \f^{y,B)\  =  const  and  the  functions  |/a(y)|  are 

shown  in  Figs.  20  and  21  for  various  incidence  angles.  As  contrasted  to  the  pressure 
wave,  the  mass-flux  disturbance  has  a  peak  near  the  upper  boundary-layer  edge.  Its 
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maximum  corresponds  to  the  incident  angle  0  ~  40° .  Important  that  this  peak  overlaps 
with  the  second-mode  peak  in  the  critical-layer  (see  Fig.  9).  Such  an  overlapping  can 
increase  the  scalar  product  <  G(x^q).F^,Bq  >  and  enhance  receptivity. 

Figure  22  shows  contours  of  the  maximum  streamwise  mass-flux,  max|/^(y)|,  in  the 

frequency-Reynolds  number  plane.  It  is  seen  that  the  mass-flux  amplitude  monotonically 
decreases  with  R  and  F . 

Fast  waves 

Characteristics  of  fast  acoustic  waves  were  calculated  at  the  same  frequency 
F  =  8.  X 10“^  and  Reynolds  number  R  =  1000 .  The  pressure  amplitude  distributions  are 
shown  in  Figs.  23  and  24.  Fast  waves  have  no  sonic  lines  and  quiet  zone  because  their 
phase  speed  is  +  a  .  Acoustic  rays  can  easily  penetrate  into  the  boundary  layer 

that  changes  the  pressure  amplitude  pattern  (compare  Figs.  18  and  23).  Because  of  this, 
the  maximum  of  IPaCy)!  is  observed  near  the  wall.  Its  extremum  corresponds  to  the 

incident  angle  0  =  50°. 

The  streamwise  mass-flux  contours,  \f^{y,&)\  =  const ,  and  the  functions,  |/a(y)|,  are 

shown  in  Figs.  25  and  26  at  various  incidence  angles.  The  mass-flux  amplitude  has  a 
peak  near  the  upper  boundary-layer  edge.  It  is  higher  and  sharper  than  in  the  slow-wave 
case  (see  Figs.  20  and  21).  Figure  27  shows  the  contours  of  max|/^(y)|  in  the  frequency- 

Reynolds  number  plane.  Similar  to  the  case  of  slow  waves  (see  Fig.  22),  the  streamwise 
mass-flux  amplitude  of  fast  waves  monotonically  decreases  with  R  and  F  . 

Foregoing  numerical  examples  indicate  that  hypersonic  boundary  layer  is  more  sensitive 
to  fast  acoustic  waves  than  to  slow  ones.  That  is  consistent  with  the  acoustic  field 
topology  predicted  by  the  short-wave  theory  (Fig.  17).  However  the  slow-wave  phase 
speed  is  closer  to  the  second-mode  phase  speed  that  is  more  preferable  for  receptivity. 
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The  last  argument  is  important  for  receptivity  on  smooth  surfaces  such  as  the  leading- 
edge  receptivity.  However,  it  is  not  critical  for  receptivity  to  sound  plus  wall  roughness, 
because  the  roughness  shape  can  be  always  tuned  to  meet  the  synchronism  condition 
(2.12). 

Small  incidence  angles 

Figures  18,  20  (for  slow  waves)  and  Figures  23,  25  (for  fast  waves)  show  that  the 
disturbance  amplitude  decreases  rapidly  as  the  incidence  angle  tends  to  zero  (phase 
speed  ^  1  ±  1  /  M ).  This  is  due  to  the  fact  that  the  reflected  acoustic  wave  is  almost  in 

counter-phase  with  the  incident  wave  and  they  cancel  each  other  in  the  boundary  layer. 
The  limit  0->O  needs  to  be  analyzed  carefully  because  the  length  scale  of  the 
amplitude  variation  in  vertical  direction  tends  to  infinity,  I*  =  (In  /  kjS^ oo . 
According  to  the  analysis  of  [2]  nonparallel  effects  can  change  the  acoustic  field 
essentially  if  the  length  ratio  is  /* /A*^  =0(A^F*'-)  (A^  =A*J^v*x*/U*  is 

nondimensional  displacement  thickness,  which  is  A^  *  14  in  the  case  considered). 
Using  Eqs.  (3.5),  (3.6)  we  obtain  that  the  local-parallel  approximation  is  not  valid,  if  the 
incident  angle  Q<0{AF^'-).  For  the  frequency  parameter  F  =  8.xl0"^  this  gives  the 
restriction  0<7°. 

The  asymptotic  analysis  [2]  indicates  that  the  boundary  layer  growth  (nonparallel  effect) 
causes  diffraction  of  incident  acoustic  waves.  As  a  result,  acoustic  field  amplitude  at  the 
upper  boundary-layer  edge  decreases  in  downstream  direction.  The  boundary-layer 
displacement  leads  to  formation  of  “quiet”  zone  near  the  plate  surface  and  far  from  the 
leading  edge. 

For  configuration  shown  in  Fig.  16,  a  major  portion  of  the  plate  surface  is  radiated  by 
acoustic  waves  of  small  incidence  angle.  This  is  typical  for  acoustic  sources  to  be  at  rest. 
Both  the  diffraction  and  the  counter-phase  interference  can  essentially  reduce  the  sound 
amplitude  in  the  boundary  layer.  This,  in  turn,  reduces  the  normal-wave  excitation  due  to 
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sound  scattering  by  a  roughness  element  located  far  downstream  from  the  leading  edge. 
However,  in  many  practical  cases,  sound  waves  are  generated  by  sources  convected 
downstream.  Their  incidence  angles  are  normally  large.  For  example,  the  measurements 
of  Laufer  [24]  indicated  that  a  turbulent  boundary  layer  on  walls  of  a  conventional  wind 
tunnel  radiates  sound  with  the  incidence  angle  larger  than  45° .  These  arguments  indicate 
that  receptivity  experiments  should  cover  a  wide  range  of  the  acoustic-wave  incidence 
angle.  Special  techniques  need  to  be  developed  to  simulate  sound  from  convective 
acoustic  sources. 

3.5  Distributed  receptivity  on  wavy  surface 

Using  Eq.  (2.20)  we  express  the  initial  amplitude  of  the  second-mode  wave  generated  by 
wavy  wall  in  the  resonance  region  as 


h2(^io’0)|  =  (3.9) 

where  |p2(^io’^)|  pressure  amplitude  on  the  wall  surface;  the  receptivity  function  is 
expressed  in  the  form 
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(3.10) 

(3.11) 


The  acoustic-wave  vector-function,  F^,  is  normalized  by  the  condition  (pressure 
amplitude  of  the  incident  wave  at  the  upper  boundary-layer  edge)=l;  i.e.  /  P*  at 

y  =  .  The  waviness  amplitude  is  measured  in  units  of  the  local  scale 

where  jCq  is  resonance-point  coordinate;  i.e.  h  =  h*  I yjvlxl  / f/*  .  The  receptivity 
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function  is  considered  as  a  function  of  frequency  parameter  F ,  Reynolds  number 
^  acoustic- wave  incidence  angle  0 . 

Because  |lm(a2)|«|Re(a2)|  for  all  parameters  under  consideration,  we  use  the 
following  approximation  of  the  resonance  condition  (2.12) 

(x„=a,(x,^,(o)-a^«Re[a,(x^^,co)-aJ.  (3.12) 

Calculations  of  the  distributed  receptivity  function  G^^,^  were  performed  for  slow  and 

fast  acoustic  waves  of  frequencies  F  =  4.xlO■^6.xlO‘^8.xlO■^  and  lO"*  in  the 
Reynolds  number  ranges  relevant  to  the  second-mode  instability  (see  Figs.  5  and  6). 

For  slow  acoustic  waves,  the  wavenumber  a„  is  expressed  by  Eq.  (3.6)  as  a  function  of 
the  incident  angle  0.  Figures  28a-31a  show  the  receptivity  function  Gp^{Q,R). 
Figures  28b-31b  represent  the  contours  G =  const .  In  all  cases,  the  receptivity 

function  is  highly  non-uniform  with  respect  to  the  Reynolds  number.  It  has  two  sharp 
maximums:  the  first  is  observed  near  the  lower  neutral  branch  and  the  second  near  the 
upper  neutral  branch.  Similar  to  the  case  of  wall-induced  disturbances  [4],  these 
maximums  correlate  with  the  discrete-spectrum  branch  points  shown  in  Rg.  8. 

The  first  maximum  is  of  the  most  interest  to  the  transition  prediction  methodology 
because  it  corresponds  to  the  initial  phase  of  instability  amplification.  The  contour  plots 
indicate  that  this  maximum  is  observed  at  relatively  high  incidence  angles,  0  «  65° ,  for 
all  frequencies  considered.  The  second  maximum  may  be  of  a  particular  interest  to 
bypass  modeling.  It  is  slightly  larger  than  the  first  maximum  and  observed  at  smaller 
incidence  angles  0  =  40°.  Receptivity  to  acoustic  waves  of  small  0  is  on  the  order  of 
magnitude  lower  than  its  maximum  value.  This  is  consistent  with  the  fact  that  such 
acoustic  waves  do  not  penetrate  into  the  boundary  layer  due  to  the  counter-phase 
interference  between  incident  and  reflected  waves  (see  Section  3.4). 
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For  fast  acoustic  waves,  the  wavenumber  is  specified  by  Eq.  (3.7).  Figures  32a-35a 
show  the  receptivity  function  and  Figures  32b-35b  illustrate  the  contours 

Gp^iQ,R)  =  const .  Receptivity  to  fast  acoustic  waves  is  higher  than  to  slow  ones.  This 

is  consistent  with  the  fact  that  the  fast-wave  amplitude  is  larger  than  the  correspondent 
slow-wave  amplitude  (see  Section  3.4).  The  first  and  second  maximums  of  the 
receptivity  function  merge  into  one  maximum  as  the  frequency  increases.  For 
F  >6.  xlO"*,  they  form  a  relatively  flat  peak  covering  the  unstable  Reynolds  number 
range.  The  receptivity  function  maximum  is  observed  at  the  incidence  angle  0  =  45° . 
This  maximum  is  a  weak  function  of  the  disturbance  frequency.  As  the  incidence  angle 
tends  to  zero,  the  receptivity  function  decreases  rapidly. 

3.6  Local  receptivity  on  hump 

Using  Eq,  (2.24)  we  can  express  the  second-mode  amplitude  generated  by  acoustic 
waves  on  a  local  hump  as 


|/>2(Xjo,0)|  —  EghG p Q),cc^)^p((x^q']^  .  (3.13) 

Here  the  local  receptivity  function  is 

GpM  =  1^)^  %4  (^10  ’0)|  >  (3-14) 


the  Fourier  component  p(a,^)  is  determined  as 

-H» 

P(ctw)=  j  g(x')exp(-ia^x')dx' ,  (3.15) 

— oo 

the  function  q  is  given  by  Eq.  (3.11).  The  resonance  wavenumber  is  approximated 
by  Eq.  (3.12).  Similar  to  the  case  of  receptivity  on  a  wavy  wall,  Gp^^  is  considered  as  a 
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function  of  the  frequency  parameter  F ,  Reynolds  number  R  =  ^Ju*x*  /v*  based  on  the 
hump  center  ,  and  the  acoustic-wave  incidence  angle  0. 


The  receptivity  functions  and  their  contours  =  const  are  shown 

in  Figs.  36-39  for  slow  acoustic  waves  and  in  Figs.  40-43  for  fast  acoustic  waves. 
Comparing  these  data  with  the  distributed  receptivity  data  shown  in  Figs.  28-35  we 
conclude  that  the  local  receptivity  function  is  much  lower  than  the  distributed  one, 
^pM  ~  .  However,  qualitative  behavior  of  G^j^^  is  similar  to  that  observed  for 

the  distributed  receptivity.  This  is  due  to  the  fact  that  the  ratio. 


^P4is  _  p~l/2 

r  “ 


2m 


{da^  /  dx^  _ 


(3.16) 


does  not  depend  on  characteristics  of  acoustic  wave  and  roughness.  In  the  cases 
considered,  the  eigenvalue  derivative  da^/dx^  is  a  regular  function  of  because  the 
spectrum  branch  points  are  shifted  to  the  complex  R  -plane  as  shown  in  Fig.  8.  For  this 
reason,  the  ratio  (3.16)  is  a  relatively  weak  function  of  the  streamwise  coordinate  Xjq  or, 
what  is  the  same,  the  Reynolds  number  R . 


Note  that  in  the  case  of  hypersonic  boundary  layer  on  a  cooled  plate  [1,  4],  the  branch 
points  may  be  close  to  the  real  axis.  Since  the  eigenvalue  derivative  tends  to  ±  ««>  in  the 
branch-point  vicinity  (see  Eq.  (2.25)),  the  local  behavior  of  the  functions  and 

(^p.ioc  will  be  different.  In  this  case,  we  should  use  the  receptivity  model  of  Section  2.3. 
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4.  Summary  discussion 

In  this  study,  we  extended  the  receptivity  model  developed  in  [4]  for  wall-induced 
perturbations  to  the  case  of  freestream  disturbances  scattered  by  surface  roughness.  We 
considered  the  second-mode  excitation  by  acoustic  waves  interacting  with  the  wall 
waviness  (distributed  roughness)  and  hump  (local  roughness  element).  This  problem 
includes  the  following  components; 

(a)  stationary  disturbances  generated  by  surface  roughness; 

(b)  acoustic  disturbance  resulted  from  the  interaction  of  incident  acoustic  wave  with  the 
boundary-layer  mean  flow  and  wall; 

(c)  unsteady  perturbation  resulted  from  the  nonlinear  interaction  of  the  roughness-induced 
disturbance  with  the  acoustic  wave. 

As  contrasted  to  the  subsonic  case,  a  2-D  acoustic  source  in  supersonic  flow  generates 
slow  waves  with  the  phase  speed  0<c^<l-l/M,  fast  waves  with  >  1  +  1/M  and 

counter-flow  waves  with  <0 .  The  latter  stagnate  near  the  Mach  line  and  form  a  local 
disturbance  with  zero  phase  speed.  In  this  report,  we  analyzed  receptivity  to  slow  and  fast 
acoustic  waves  only. 

It  is  shown  that  the  receptivity  mechanism  depends  on  the  discrete  spectrum  topology.  If 
the  spectrum  is  simple  (i.e.  its  branch  points  are  essentially  complex  and  their  influence 
domains  do  not  overlap  with  the  real  parameter  space),  then  the  receptivity  mechanism 
reveals  the  following  features; 

•  For  distributed  roughness,  a  strong  excitation  occurs  in  a  local  region  of  the  length 
/*  =  ,  where  Component  (c)  is  in  resonance  with  the  second-mode 

wave.  In  this  region,  the  normal-wave  amplitude  is  proportional  to  (ZT  / . 

As  the  Reynolds  number  increases,  the  relative  length  of  the  resonance  region,  V 16* , 
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increases.  As  /?  ->  oo ,  the  normal  wave  amplitude  tends  to  infinity.  In  this  case,  the 
upper  limit  of  the  disturbance  amplitude  is  determined  by  nonlinear  effects  of  the 
higher-order  approximation. 

•  For  local  roughness,  the  normal-wave  amplitude  is  finite  even  in  the  parallel  flow  limit. 
This  is  due  to  the  fact  that  only  a  small  portion  of  the  roughness  spectrum, 
{cc„- a„Q)  =  OiR~^'^),  is  involved  into  the  resonance  mechanism.  As  J?— >oo,  the 

resonant  excitation  increases  proportionally  to  R}'^,  whereas  the  active  spectrum 
region  is  narrowed  down  proportionally  to  R~''^ .  Because  both  trends  compensate 
each  other,  the  normal-wave  amplitude  tends  to  a  finite  value. 

If  the  branch-point  influence  domain  overlaps  with  the  real  parameter  space,  then  both 
propagation  and  excitation  of  normal  waves  are  singular.  The  receptivity  function  tends  to 
mfinity  as  the  resonance  point  approaches  the  spectrum  branch  point.  Asymptotic  analysis 
of  this  singularity  leads  us  to  the  following  conclusions; 

•  The  second-mode  excitation  on  a  wavy  wall  (distributed  roughness)  is  not  singular  at 
the  branch  point.  The  receptivity  function  peak  is  totally  due  to  a  singular  behavior  of 
the  normal  wave  itself  If  one  measures  the  disturbance  amplitude  (2  at  a  fixed  point 

downstream  fi’om  the  branch  point  for  various  loci  of  the  resonance  point  x,o , 
then  he  will  not  observe  essential  changes  in  the  distribution  2(x,o)  near  the  point 
^10  ~  ^16  • 

•  The  second-mode  excitation  over  a  hump  (local  roughness)  strongly  depends  on  a 
distance  between  the  hump  locus  x,o  and  the  branch  point.  As  Xjj ,  the  normal- 

wave  amplitude  increases  proportionally  to  |x,o  -Xj^p'^'^and  attains  its  maximum  value 
of  the  order  of  at  x,o  =  x,j, .  If  one  measures  the  disturbance  amplitude  0  at 

a  fixed  point  downstream  from  the  branch  point  for  various  loci  of  the  forcing  element, 
then  he  will  observe  a  local  peak  in  the  distribution  2(^io)  at  the  point  x,o  =  x,j, . 
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The  analytical  results  were  build  in  a  computational  module,  which  was  coupled  with  our 
stability  solver.  Detailed  calculations  were  performed  for  two-dimensional  disturbances 
interacting  with  the  flat-plate  boundary  layer  at  the  Mach  number  M  -6,  wall 
temperature  J  *=  80°  F  and  stagnation  temperature  =  390°  F.  This  case  is  considered 

as  a  good  candidate  for  receptivity  experiments  in  the  Mach  6  quiet  wind  tunnel.  The 
spectrum  calculations  showed  that  the  branch  points  are  essentially  complex.  That 
motivated  us  to  perform  receptivity  calculations  under  the  assumption  that  the  spectrum  is 
simple.  Numerical  studies  lead  us  to  the  following  conclusions: 

1 .  In  the  boundary  layer,  the  pressure  amplitude  of  fast  acoustic  waves  is  larger  than  that 
of  slow  waves  with  the  same  frequency  and  incidence  angle.  The  amplitude  maximum 
is  observed  at  the  incidence  angle  0  «  40°  for  slow  waves  and  0  «  45°  for  fast  waves. 
Due  to  counter-phase  interference  between  the  incident  and  reflected  waves,  the 
disturbance  amplitude  decreases  rapidly  as  the  incidence  angle  tends  to  zero.  This 
trend  is  opposite  to  the  subsonic  case  when  the  incident  and  reflected  waves  are  in 
phase. 

2.  For  receptivity  on  wavy-wall,  the  receptmty  function  is  highly  non-uniform  with 
respect  to  the  Reynolds  number. 

•  In  the  case  of  slow  acoustic  waves,  it  has  two  maximums:  the  first  is  close  to  the 
lower  neutral  branch  and  the  second  to  the  upper  neutral  branch.  These  maximums 
correlate  with  the  spectrum  branch  points.  The  boundary  layer  is  the  most 
receptive  to  acoustic  waves  of  relatively  high  incidence  angles;  0  «  65°  near  the 
first  maximum  and  0  »  40°  near  the  second  maximum. 

•  In  the  case  of  fast  acoustic  waves,  the  receptivity  maximums  merge  into  one 
relatively  flat  maximum  covering  the  instability  region.  The  boundary  layer  is  most 
receptive  to  acoustic  waves  of  the  incidence  angle  0  «  45° . 
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•  Receptivity  to  both  slow  and  fast  acoustic  waves  of  small  incidence  angles  is  on 
the  order  of  magnitude  lower  than  its  maximum  level. 

3.  For  receptivity  on  local  hump,  the  receptivity  function  has  the  same  qualitative 
behavior  as  in  the  wavy-wall  case.  However  its  quantity  is  essentially  smaller: 


This  study  does  not  cover  all  aspects  of  the  hypersonic  boundary-layer  receptivity  to 

sound  scattered  by  roughness.  The  following  issues  \\ill  need  to  be  addressed  in  future 

effort: 

•  For  the  boundary  layers  on  cooled  surfaces,  the  spectrum  branch  points  may  be  very 
close  to  the  real  parameter  space.  Careful  studies  of  the  receptivity  mechanism  in  the 
branch  point  vicinity  should  be  performed. 

•  Calculations  have  been  conducted  for  two-dimensional  cases  only.  Since  in  the 
majority  of  practical  cases  (including  wind-tuimel  experiments),  acoustic  sources  and 
roughness  elements  are  three-dimensional,  analysis  of  3-D  receptivity  will  be  useful. 

•  Receptivity  to  acoustic  disturbances  concentrated  near  the  line  of  intersection  between 
the  Mach  cone  and  body  surface  may  play  important  role  in  the  instability  excitation. 
This  problem  seems  to  be  relevant  to  the  boundary-layer  interaction  with  oscillating 
shocklets. 

•  Parametric  studies  of  hypersonic  receptivity  to  entropy  and  vorticity  disturbances 
scattered  by  surface  irregularities  are  also  required. 

•  It  is  feasible  to  compare  the  results  of  this  study  with  our  previous  results  relevant  to 
the  leading-edge  receptivity  [1],  vorticity/entropy  swallowing  [2],  and  receptivity  to 
wall-induced  perturbations.  This  will  help  to  identify  the  most  probable  scenarios  of 
the  initial  transition  phase  for  typical  practical  cases. 
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Fig.  1  Flow  scheme. 


free  stream 


Fig.  2  Flow  scheme  for  resonant  excitation. 
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Fig.  3  Mean-flow  velocity  profiles  U(y),  and  dU  Idy . 
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Fig.  4  Mean-flow  temperature  profiles  T{y),  and  dT/dy . 
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Fig.  6  Second-mode  growth  rate  -  Im(a2 )  as  a  function  of  Reynolds  number  at  various 
frequencies. 


49 


Fig.  8  Branch  points  in  complex  R  -plane  at  various  frequencies. 
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Fig.  9  Second-mode  eigenfunctions  for  pressure  streamwise  mass  flux 

near  lower  neutral  curve:  R  =  1000 ,  F  =  8.  x  10"^ . 
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Fig.  10  Contours  Re{p„(y,a„))  =  const  for  wall-induced  disturbances;  /?  =  1000. 


Fig.  11  Distributions  Re(/7„(y))  for  wall-induced  perturbation  at  various  wavenumbers 
a  =  a^  ,  R  =  1000. 
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Fig.  12  Contours  Re(/„(y,a^))  =  con^t  for  wall-induced  disturbace;  /f  =  1000 
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Fig.  13  Mass-flux  distributions  Re(/„(y))  for  wall-induced  perturbation  at  vanous 
wavenumbers  a  =  a,/,  R  =  1000 . 
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Fig.  14  Pressure  amplitude  contours  |p„,(>',a)|  =  const ;  R  =  1000 . 


Fig.  15  Mass-flux  amplitude  contours  |/„(y,a)|  =  const ;  R  =  1000 . 
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Fig.  16  Acoustic  waves  from  cylindrical  source. 


Fig.  17  Scheme  of  the  acoustic  wave-boundary  layer  interaction. 


Fig.  21  Distributions  \faiy,&)\  at  various  incidence  angles  for  slow  waves  /?  =  1000 
F  =  S.xlO-\ 
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Fig.  23  Contours  of  pressure  amplitude  \p^{y,Qi)\  =  const  for  fast  waves;  /?  =  1000, 
F  =  8.xl0-\ 


Fig.  24  Distributions  |Pa(y’€))|  various  incidence  angles  for  fast  waves,  R  =  1000, 
F  =  8.xl0-^ 
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Fig.  27  Contours  max|/„(>',  R,  F)|  =  const  for  fast  waves  at  0  =  45° . 
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Fig.  28a  Distributed  receptivity  function  for  slow  acoustic  waves;  F  =  4.  x  10 


Fig.  28b  Contours  G  ^,^(0,  R)  =  const  for  slow  acoustic  waves;  F  =  4.  x  10 


Fig.  29a  Distributed  receptivity  function  for  slow  acoustic  waves;  F  =  6.  x  10 
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Fig.  29b  Contours  ^,.^(0,  R)  =  const  for  slow  acoustic  waves;  F  =  6.  x  10  ^ 


Fig.  30a  Distributed  receptivity  function  for  slow  acoustic  waves;  F  =  8.  x  10 
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Fig.  30b  Contours  G.  ....(0,  R)  =  const  for  slow  acoustic  waves;  F  = 
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Fig.  37a  Local  receptivity  function  for  slow  acoustic  waves;  F  =  6.  x  10  ^ 
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Fig.  37b  Contours  R)  =  const  for  slow  acoustic  waves;  F  =  6.x  10  ^ 
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Fig.  38a  Local  receptivity  function  for  slow  acoustic  waves;  F  =  8.  x  10  ^ . 
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Fig.  38b  Contours  (0,  R)  =  const  for  slow  acoustic  waves;  F  =  8.  x  10 
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Fig.  39b  Contours  Gpj^^{Q,R)  =  const  for  slow  acoustic  waves;  F  =  10”^ . 
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Fig.  40a  Local  receptivity  function  for  fast  acoustic  waves;  F  =  4.  x  10  * 
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Fig.  40b  Contours  R)  =  const  for  fast  acoustic  waves;  F  =  4.  x  10 


Fig.  41a  Lx)cal  receptivity  function  for  fast  acoustic  waves;  F  =  6.xl0  ^ 
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Fig.  41b  Contours  R)  =  const  for  fast  acoustic  waves;  F  =  6.x  10  ^ 


Fig.  43a  Local  receptivity  function  for  fast  acoustic  waves;  F  =  10  . 
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Fig.  43b  Contours  Gpj^(@,R)  =  const  for  fast  acoustic  waves;  F  =  10"^ . 


Appendix 

Non-zero  elements  of  Matrix  Hq  =  {aj,),  j  =  ,  /  =  1,. 

^12  “  ^56  ~  ^78  ~  ^ 

a,,  =a'  +  P'  +/?-^ 

^  nT 


Cl^'i  — 


DJI 


DT 


DH^^DU 


a,3  =  -ia(m  -t- 1) - ioc  — —  +  R 


HT 


a,,  =  R—-(m  +  l)yM~(xQ 
1^ 

,  ,3^  D(ii'DU) 

a,5  =a(m  +  l) - - - 


^26  ~ 


H'DU 
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=-ia 


DT 


«33=- 


=-iyM'Q. 
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«35 

^37  ~  ~^P 


«4i  =-«Za 


r.  51+ 25a 

T  /i 


«42  = 
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^43  ~  ^ 


iQ 


D-T  DpiDT 


.(X^-  _R^-  ^R:r:L+r—+r 

^  llT  T  llT 


^44  =-iXryM- 


ocDU  +  PDW  +  Q 

f  DT  DjU  Y 

«45 


«46  =iXr— 


{aDU  +  pDwi^j+^ 
Q 


+  r 


QZ)jU 


^47  = 


«48  =-ixP 


DT^.Dij. 

r - +  2—^ 

T  M 


..8  in  Eq.  (2.9): 
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where  k  =  0  corresponds  to  Stokes  hypothesis; 

U (x,,  y)  and  y)  is  x  and  z-component  of  the  laminar  mean  flow  on  a  swept  wing; 
r(x, ,  y)  is  mean-flow  temperature; 

jU  =  fi(T)  is  nondimensional  dynamic  viscosity,  D/t  =  n'DT . 
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Non-zero  elements  of  the  nonlinear  interaction  matrix 
G(y,x^,q;a„,P„,a^,p^,co)  =  (gji),  j  =  in Eq.  (2.8a): 


^21  +CcJqi+Pa^j+TQ^p] 


g,,=—(q,+TDUp) 
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832  =Y  ^a.^3 - ^^5+^6 


834  =iyM-{-a„q,  -p„q^  +TCl^p) 

835  ~  ■*■  ^0^7 )  “  ^^awP  j,2  ^3 


^36  ““  j, 

837  =  -iPwTp 

■CCa 

841=1  j  •■<13 

1  r  DT 

843  =  2—^3  +  i(P^  -  Pjq,  -I- iT[(a^  -aJU  +  - PJW  - (o]p 

Su  =  -oc.w+(.^.-^.)w-  «>]«, 

*.J  =  -  ^[(a,  -  a,  W  +  ( A  -  )W'  -  ®]9j 

g,T=j<.P.-Pj4, 

g„=ia,Pr-^[q,-(y-l)M‘TqJ 
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